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Abstract 

The decay widths of top quark to S-wave be and bb bound states are evaluated at the next- 
to-leading(NLO) accuracy in strong interaction. Numerical calculation shows that the NLO cor- 
rections to these processes are remarkable. The quantum chromodynamics(QCD) renormalization 
scale dependence of the results is obviously depressed, and hence the uncertainties lying in the 
leading order calculation are reduced. 
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I. INTRODUCTION 



Since predicted by the Standard Model (SM) , top quark has become an important 
role in high energy physics due to its large mass, which is close to the electroweak symmetry 
breaking scale [4j. A great deal of researches focusing on top quark physics have been 
performed after its discovery in 1995 in the Fermilab |5|. On the experiment aspect, with 
the running of the Tevatron and forthcoming LHC, the lack of adequate events will not be an 
obstacle for the top quark physics study According to Ref . [6| , at the LHC 10 7 ~ 10 8 ti pairs 
can be obtained per year, so this enables people to measure various top quark decay channels. 
Meanwhile, the copious production of the top quarks supplies also a great number of bottom 
quark mesons since the dominant top quark decay channel is t — > b + W + . Therefore, 
the bottom quark meson production in top quark decays may stand as an important and 
independent means for the study of heavy meson physics and the test of perturbative QCD 
(pQCD). 

As the known heaviest mesons, bottomonia and B c (be or be) possess particular meaning 
in the study of heavy flavor physics. The LHCb as a detector specifically for the heavy 
flavor study at the LHC will supply copious B c and T data for this aim. Theoretically, the 
direct hadroproduction of B c and T was studied in the literature [t 9 ] . In addition to the 
"direct" production, "indirect" process as in top quark decays may stand as an independent 
and important source for B c and T production. Since the top quark's lifetime is too short 
to form a bound state [10], the B c and T production involved scheme in top quark decays 
is less affected by the non-pertubative effects than in other processes. In Ref. ll), the top 
quark decays into T and B* at the Born level was evaluated. Recently, the S- and P-wave 
B c meson productions in top quark decays were fully evaluated, including the color-octet 
contributions, at the leading order accuracy of QCD by Chang et al. [12 ]. 

Considering the importance of investigating B c and T in the study of perturbative Quan- 
tum Chromodynamics (pQCD) and potential model, it is reasonable and interesting to eval- 
uate the production rates of these mesons in top quark decays at the next-to-leading order 
(NLO) accuracy of pQCD. At the bottom quark and charm quark mass scales the strong 
coupling is not very small, therefore the higher order corrections are usually large. On the 
other hand, in the processes of top quark decays into B C (Y), the t —> bc(b) + c(b) + W + , 



there exist large scale uncertainties in the tree level calculation [13|. The NLO corrections 
should in principle minimize it and give a more precise prediction. To calculate the B c and 
T production rates in top quark decays at the NLO accuracy are the aims of this work. In 
our calculation, both of the S-wave spin-singlet and -triplet states are taken into account, 
i.e., B*, B C) T and rjb- To deal with the non-perturbative effects, the non-relativistic QCD 
(NRQCD) 14| effective theory is employed. The calculation will be performed at the NLO 
in pQCD expansion, but at leading order in relativistic expansion, that is in the expansion 
of v, the relativistic velocity of heavy quarks inside bound states. 

The paper is organized as follows: after the Introduction, in section II we explain the 
calculation of leading order decay width. In section III, virtual and real QCD corrections 
to Born level result are evaluated. In section IV, the numerical calculation for concerned 
processes at NLO accuracy of pQCD is performed, and the scale dependence of the results 
is shown. The last section is remained for a brief summary and conclusions. 



II. CALCULATION OF THE BORN LEVEL DECAY WIDTH 

At the leading order in a s , there are two Feynman Diagrams for each meson production, 
which are shown in Figure 1. For the convenience of analytical calculation, taking B c as an 
example, the momentum of each particle is assigned as: p\ = p t , p% = Pb, Pa = Pc, P5 = p c , 
Pq — Vw+i Po — P3 + Pi, P3 — —P4- For bottomonium, the only difference is that 04 and 

Trie 

p 5 represent the momenta of anti-bottom quark and bottom quark which are produced in 
gluon splitting. 

Of the B c and B* production in top quark decays, i.e. 

tip,) -+ B c /B* c (p ) + c(p 5 ) + W + (p 6 ) , (1) 
we employ the following commonly used projection operators for quarks hadronization: 



and 





(a) (b) 

FIG. 1: The leading order Feynman diagrams for B c and B* production in top quark decays. 

Here, £#* is the polarization vector of B* with po ■ e = 0, l c stands for the unit color 
matrix, and N c = 3 for QCD. The nonperturbative parameters ip§ c (0) and ^g*(0) are the 
Schrodinger wave functions at the origin of be bound states, and in the non-relativistic limit 
V'BcCO) = m our calculation, the non-relativistic relation m^ c = mg, = mj + m c is 

also adopted. 

The LO amplitudes for B c production can then be readily obtained with above prepara- 
tions. They are: 

-^a = 7== — u{p B )j^5{p0 + msjj/i-f ; 2^ j \2 u \P l ) . 4 



and 



7ra a gipB c ( ) v tb c F S 3,k-/ . , _ ^(pe)(l - 7s) (^3+ jfe + mt) 



M 6 = 7== -uipsh^si^o + ™ Bc )— — - — — • (5) 



mr 



Here, j, are color indices, Cp = 4/3 belongs to the SU(3) color structure. For £>* 
production, the amplitudes can be obtained by simply substituting 275(^0 + m B c ) with 
4b* (A + m b + rn c ) in above expressions. 

The Born amplitude of the processes shown in Fig.l is then M.Bom = -Ma + -Mb, an d 
subsequently, the decay width at leading order reads: 

dTB ° rn = 2^~2 iV l-^Born| 2 dPS 3 • (6) 

Here, represents the sum over polarizations and colors of the initial and final particles, \ 
and come from spin and color average of initial t quark, dPS3 stands for the integrants 



of three-body phase space, whose concrete form is 

dPS 3 = J^-i^d Sl d5 2 , (7) 

where Si = (P0+P5) 2 = (pi— Po) 2 and s 2 = (j>5+Pq) 2 = (pi—Po) 2 are Mandelstam variables. 
The upper and lower bounds of the above integration are 



max 
£>1 



ijf[m 2 , s 2 , m|J • f[s 2 , m 2 , rri^] + [mf - s 2 - (m b + m c ) 2 ] (s 2 + m 2 - m^) 
+m 2 s +m 2 c , (8) 
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f[m 2 ,s 2} m 2 B J ■ f[s 2l m 2 cl m 2 w \ - [m 2 - s 2 - (m b + m c ) 2 ](s 2 + m 2 



in 



w) 



2s 2 

,2 , 2 



+m 2 Bo + mi (9) 
and 

sT x = [m t - (m b + m c )} 2 , = K + m w ) 2 (10) 

with 

f[x,V,z] = -y - z) 2 - Ayz . (11) 

III. THE NEXT-TO-LEADING ORDER CORRECTIONS 

At the next-to-leading order, the top quark decays to B c and T include the virtual 
and real QCD corrections to the leading order process, as shown in FigsfSE With virtual 
corrections, the decay widths at the NLO can be formulated as 

dT Virtual = JL-_ J] 2Re(M* Born M Virtual )dPS 3 . (12) 

The ultraviolet (UV) and infrared(IR) divergences usually exist in virtual corrections. We 
use the dimensional regularization scheme to regularize the UV and IR divergences, similar 



as performed in Ref. 



15], and the Coulomb divergence is regularized by the relative velocity 



v. In dimensional regularization, 75 is difficult to deal with. In this calculation, we adopt 
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FIG. 2: The self-energy diagrams in virtual corrections. 



the Naive scheme, that is, 75 anticommutates with each 7^ matrix in d-dimension space- 
time, {75, 7 M } = 0. The UV divergences exist merely in self-energy and triangle diagrams, 
which can be renormalized by counter terms. The renormalization constants include Z 2 , 
Z 3 , Z m , and Z g , corresponding to quark field, gluon field, quark mass, and strong coupling 
constant a s , respectively. Here, in our calculation the Z g is defined in the modified-minimal- 
subtraction (MS) scheme, while for the other three the on-shell (OS) scheme is adopted, 
which tells 



Here, /3q = (11/3)C^ — (4/3)TfUf is the one-loop coefficient of the QCD beta function; 
rif = 5 is the number of active quarks in our calculation; Ca = 3 and Tp — 1/2 attribute to 
the SU(3) group; fi is the renormalization scale. 

In virtual corrections, IR divergences remain in the triangle and box diagrams. Of 
all the triangle diagrams, only two have IR divergences, which are denoted as TriangleN7 
and TriangleN9 in Figj3j Of the diagrams in FigJU BoxN3 has no IR singularity, while 
BoxN4 and PentagonN9 have Coulomb singularities and PentagonN9 possesses ordinary 





- — — 7£ + ln47r + 0(e) 
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(13) 
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FIG. 3: The triangle diagrams in virtual corrections. 



IR singularity as well. The remaining diagrams all have IR singularities, while the combi- 
nations BoxN2 + BoxN6, BoxNl + PentagonN8 + TriangleN9, BoxN5 + TriangleN7 are IR 
finite. The Coulomb singularities belonging to BoxN4 and PentagonN9 can be regularized 
by the relative velocity v. After regularization procedure, the \ term will be canceled out 
by the counter terms of external quarks which form the B c or T, while the ^ term will be 
mapped onto the wave functions of the concerned heavy mesons. The remaining IR singu- 
larities in BoxN7 and BoxN9 are canceled by the corresponding parts in real corrections. In 
the end, the IR and Coulomb divergences in virtual corrections can be expressed as 



I R,Coulomb 



dr 

virtual 



dr# orri — — 

671 



V 



1 Pt ■ PcX s In x s 1 



(14) 



with p t = pi, p c = p 5 and x s = 



e m c m t (l — x 2 s ) e 

1- \J 'l-2m c mt / '(m c m t -pfPc) TT , n . i 1 • r 

— v , . Here, m this work - m ract 

l+y/l-2m, c m t /(m c m t -pfPc) 



represents \ — ^ E + ln(47r/i 2 ). 

Of our concerned processes, there are 12 different diagrams in real correction, as shown 



BoxN1 



Box N2 



BoxN3 






Box N4 



BoxN5 



Box N6 






BoxN7 



Pentagon N8 



Pentagon N9 



FIG. 4: The box and pentagon diagrams in virtual corrections. 



in FigEJ Among them, RealN2, RealN3, RealN8, and RealN9 are IR-finite, meanwhile the 
combinations of RealNl + RealN5 and RealNIO + RealNll exibit no IR singularities as well, 
due to the reasons of gluon connecting to the b or c quark of final B c or T. The remaining 
diagrams, RealN4, RealN6, RealN7, and RealN12 are not IR singularity free. To regularize 
the IR divergence, we enforce a cut on the gluon momentum, the p^. The gluon with energy 
Py < 5 is considered to be soft, while p® > 5 is thought to be hard. The 5 is a small quantity 
with energy- momentum unit. In this case, the IR term of the decay width can then be 
written as: 



where dPS4 is the four-body phase space integrants for real correction. Under the condition 
of p® < 5, in the Eikonal approximation we obtain 




(15) 
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FIG. 5: The real correction Feynman diagrams that contribute to the production of B c or B*. 
In the small 5 limit, the IR divergent terms in real correction can therefore be expressed as 
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- - Log(«5 2 
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p t ■ p c x s In x s 
m c m t (l — xl) 



+ finite terms 



(17) 



Here, the Log(<5 2 ) involved terms will be canceled out by the ^-dependent terms in the hard 
sector of real corrections. Referring to the Eq. (ll4j) . it is obvious that the IR divergent 
terms in real and virtual corrections cancel with each other. In case of hard gluons in real 
correction, the decay width reads 



2m t 2 N t 

In this case the phase space dPS 4 \hard can be written as 
2 (sy + m 2 c — m^) 2 — ^sym 



Real] 2 dPS4 \hard 



(18) 
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P7- rv+ rP7+ py+ 

dp 7 ° dy+ dp 7 ° 



PT- 



dy 



(19) 



with 



Po_ = m b + m c 



(20) 



s + m 2 b - + 2m h m c - 2m w ■ m c 
Po + ^ , (21) 

o _ s + ml - m 2 w + 2m b m c - 2m w ■ m c - 2 v /ip ° 
'"' 2 v ^-2p ° + 2 v ^M 

o _ s + m 2 b - mly + 2m b m c - 2m w ; rn c - 2^/sp Q ° 

P7+ ~ 2^-2p^-2y/W\ ' 

V- = - Po° - Pr ) 2 - l^l 2 - (P7°) 2 - 2|^b° 7 ] , (24) 

y+ - - Po ° - P7°) 2 - l^l 2 - (P7°) 2 + 2|^|p 7 °] , (25) 
where y is a dimensionless parameter defined as y — (pi — po — pi) 2 j s with y/s = m t , and 



1^1 = V^ ) 2 - m k- ( 26 ) 

The sum of the soft and hard sectors gives the total contribution of real corrections, i.e., 

■p ~pIR I j^hard 

Real — 1 Real "r 1 i?eaZ ■ 

With the real and virtual corrections, we then obtain the total decay width of t quark 
to B c and T at the NLO accuracy of QCD 

I" ^ total I" Bom ~\~ ^Virtual ~\~ ^ Real + 0(af) . (27) 

In above expression, the decay width is UV and IR finite. In our calculation the FeynArts flf] ] 



was used to generate the 



17J, and the LoopTools 



eynman diagrams, the amplitudes were generated by the FeynCalc 



18[ was employed to calculate the Passarino-Veltman integrations. 



The numerical integrations of the phase space were performed by the MATHEMATICA. 



IV. NUMERICAL RESULTS 

To complete the numerical calculation, the following ordinarily accepted input parame- 
ters are taken into account: 

m b = 4.9 GeV, m c = 1.5 GeV, m t = 17 '4 GeV, m w = 80 GeV , (28) 
ijj s (0) = Vs*(0) = ^0- = 0.3616 GeV 3/2 , (29) 

V47T 



^O(O) = rf°(0) = ^B- = 0.6812 GeV 3 / 2 , (30) 

v4tt 



tP» lo (0) = ^ LO (0) = = S (°) = 0.8277 GeV 3 / 2 , (31) 

T V ; Vb V ; V4^ V47T - lQC F a s ' V ; 



in 



TABLE I: The decay widths of the processes t -> B* + W + + c, t -)• B c + W + + c, i -)• T + W + + b 
and t — > rji, + + 6 at the tree level and with the NLO QCD corrections are presented in two 
renormalization scale \i limits, those are 2m c and mt for the first two processes and 2m b and mt 
for the other two. 





t -> 5* + w+ + c 


B c + W+ + c 


t -> T + + b 


t^r] b + W + + b 


A* 




mt 


2m c 


mt 


2rrib 


mt 


2m\j 


m t 




0.793MeV 


0.151MeV 


0.572MeV 


0.109MeV 


26.8keV 


9.54keV 


27.1keV 


9.67keV 




0.619MeV 


0.307MeV 


0.514MeV 


0.227MeV 


52.3keV 


28.2keV 


34.3keV 


24.5keV 



V tb = 1.0, Gp = 1.1660 x GeV -2 . (32) 

Here, Vtb is the Cabibbo-Kobayashi-Maskawa(CKM) matrix element and Gf is weak inter- 
action Fermi constant. 

In above numerical calculation inputs, the radial wave function at the origin for S- 



wave B*(B C ) system is estimated by potential model 19j, while the corresponding T(^) 
nonperturbative parameter is determined from its electronic decay rate [8]. One loop result 
of strong coupling constant is taken into account, i.e. 

a s(^) = --, — 5 — w — ; ,fi x • ( 33 ) 



;il - fn^LogC^ 



QCD 



With the above preparation, one can readily obtain the decay widths of top quark to be 
and bb mesons, as listed in Table [H To see the scale dependence of the LO and NLO results, 
the ratios T(fi)/T(2m c ) for be system and r(/i)/r(2m&) for bb system are showed in Figures 
[6]and[7l respectively. Calculation tells that after including the NLO corrections, the energy 
scale dependence of the results is reduced, as expected. 



V. SUMMARY AND CONCLUSIONS 

In this work we have calculated the decay widths of top quark to S-wave be and bb bound 
states at the NLO accuracy of perturbative QCD. Considering that there will be copious 
ti data in the near future at the LHC, our results are helpful to the study of the indirect 
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FIG. 6: The ratio T{fj)/T{2m c ) versus renormalization scale fi in t quark decays. The left diagram 
for the be spin-singlet state B c and the right diagram for the spin-triplet state B*. 
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FIG. 7: The ratio r(/z)/r(2m&) versus renormalization scale flint quark decays. The left diagram 
for the bb spin-singlet state 7/& and the right diagram for the spin-triplet state T. 

production of these states. They may be also useful to the future study of NLO heavy quark 
to be and bb bound states fragmentation functions. 

Numerical results indicate that the NLO corrections greatly enhance the LO results for 
bb system, while slightly decrease the be states production widths. The main reason for 
this difference is that the NLO wave function for bottomonium is much larger than that of 
LO one, while for the calculation of B c meson, the same wave function given by potential 
model is used. Although from Table I, superficially the number of indirectly produced B c 
overshoots that of T, experimentally to detect the latter is much easier than the former. 

12 



Since top quark dominantly decays into b and W + final state with a width of 1.5 GeV or 
so, numerical results remind us that the T indirect production from top quark decay is 
detectable, while it is hard to pin down the B c states by this way. 

The numerical calculation also shows that the next-to-leading order QCD corrections 
to processes t — » bc(bb) + W + + c(b) decrease the energy scale dependence of the decay 
widths as expected, and hence the uncertainties in theoretical estimation. Future precise 
experiment on the concerned processes may provide a test on the theoretical framework for 
heavy quarkonium production and the reliability of perturbative calculation for them. 
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